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Summary 

An analytical method has been developed to determine gear- 
tooth surface coordinates of face-milled spiral bevel gears. The 
method uses the basic gear design parameters in conjunction 
with the kinematical aspects of spiral bevel gear manufacturing 
machinery. A computer program entitled “SURFACE” was 
developed to calculate the surface coordinates and provide 
three-dimensional model data that can be used for finite 
element analysis. Development of the modeling method and 
an example case are presented in this report. This method of 
analysis could also be applied in gear inspection and near-net- 
shape gear forging die design. 

Introduction 

Spiral bevel gears are currently used in all helicopter power 
transmission systems. This type of gear is required to turn the 
corner from a horizontal engine to the vertical rotor shaft. 
These gears carry large loads and operate at high rotational 
speeds. Recent research has focused on understanding many 
aspects of spiral bevel gear operation, including gear geometry 
(refs. 1 to 12), gear dynamics (refs. 13 to 15), lubrication 
(ref. 16), stress analysis and measurement (refs. 17 to 21), 
misalignment (refs. 22 and 23), and coordinate measurements 
(refs. 24 and 25), as well as other areas. 

Research in gear geometry has concentrated on 
understanding the meshing action of spiral bevel gears (refs. 8 
to 11). This meshing action often results in much vibration 
and noise due to an inherent lack of conjugation. Vibration 
studies (ref. 26) have shown that in the frequency spectrum 
of an entire helicopter transmission, the highest response can 
be that from the spiral bevel gear mesh. Therefore if noise 
reduction techniques are to be implemented effectively, the 
meshing action of spiral bevel gears must be understood. 

Also, investigators (refs. 18 and 19) have found that typical 
design stress indices for spiral bevel gears can be significantly 
different from those measured experimentally. In addition to 
making the design process one of trial and error (forcing one 
to rely on past experience), this inconsistency makes 
extrapolating over a wide range of sizes difficult, and an overly 
conservative design can result. 

Research has been ongoing in an attempt to predict stresses 
(i.e., bending and contact) by using the finite element method. 

A great deal of work (refs. 27 to 30) has gone into finite 


element modeling of parallel axis gears to determine the stress 
field. Loads are typically applied at the point of highest single 
tooth contact, and then the stress in the fillet region is 
examined. Computer programs that perform this type of 
analysis are usually two dimensional in nature and have 
computer storage requirements that are small enough for 
personal computers. These attributes make them very popular 
and attractive to designers. However, a limited number of 
researchers (refs. 16 and 21) have investigated finite element 
analysis of spiral bevel gears. 

Parallel axis components (involute tooth geometry) have 
closed-form solutions that determine surface coordinates. 
These coordinates can be used as input to finite element 
methods and other analysis tools. Spiral bevel gears, on the 
other hand, do not have a closed-form solution to describe their 
surface coordinates. Coordinate locations must be solved 
numerically. This process is accomplished by modeling the 
kinematics of the cutting or grinding machinery and the 
geometry of the basic gear design. 

The objective of the research reported herein was to develop 
a method for calculating spiral bevel gear-tooth surface 
coordinates and a three-dimensional model for finite element 
analysis. Accomplishment of this task required a basic 
understanding of the gear manufacturing process, which is 
described herein by use of differential geometry techniques 
(ref. 1). Both the manufacturing machine settings and the basic 
gear design data were used in a numerical analysis procedure 
that yielded the tooth surface coordinates. After the tooth 
surfaces (drive and coast sides) were described, a three- 
dimensional model for the tooth was assembled. A computer 
program, SURFACE, was developed to automate the 
calculation of the tooth surface coordinates, and hence, the 
coordinates for the gear-tooth three-dimensional finite element 
model. The development of the analytical model is explained, 
and an example of the finite element method is presented. 

Determination of Tooth Surface 
Coordinates 

The spiral gear machining process described in this paper is 
that of the face-milled type. Spiral bevel gears manufactured in 
this way are used extensively in aerospace power transmissions 
(i.e., helicopter main/tail rotor transmissions) to transmit power 
between horizontal gas turbine engines and the vertical rotor 
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shaft. Because spiral bevel gears can accommodate various shaft 
orientations, they allow greater freedom for overall aircraft 
layout. 

In the following sections the method of determining gear-tooth 
surface coordinates will be described. The manufacturing process 
must first be understood and then analytically described. 
Equations must be developed that relate machine and workpiece 
motions and settings with the basic gear design data. The 
simultaneous solution of these equations must be done 
numerically since no closed-form solution exists. A description 
of this procedure follows. 

Gear Manufacture 

Spiral bevel gears are manufactured on a machine like the one 
shown in figure 1 . This machine cuts away the material between 
the concave and convex tooth surfaces of adjacent teeth 
simultaneously. The machining process is better illustrated in 


figure 2. The head cutter (holding the cutting blades or the 
grinding wheel) rotates about its own axis at the proper cutting 
speed, independent of the cradle or workpiece rotation. The head 
cutter is connected to the cradle through an eccentric that allows 
adjustment of the axial distance between the cutter center and 
cradle (machine) center, and adjustment of the angular position 
between the two axes to provide the desired mean spiral angle. 
The cradle and workpiece are connected through a system of 
gears and shafts, which controls the ratio of rotational motion 
between the two (ratio of roll). For cutting, the ratio is constant, 
but for grinding, it is a variable. 

Computer numerical controlled (CNC) versions of the cutting 
and grinding manufacturing processes are currently being 
developed. The basic kinematics, however, are still maintained 
for the generation process; this is accomplished by the CNC 
machinery duplicating the generating motion through point-to- 
point control of the machining surface and location of the 
workpiece. 



Figure 1 .—Machine used to generate spiral bevel gear-tooth surface. 


Cradle axis , , /— Cutter axis 




Figure 2.— Orientation of workpiece to generation machinery. 

Coordinate Transformations 

The surface of a generated gear is an envelope to the family 
of surfaces of the head cutter. In simple terms this means that 
the points on the generated tooth surface are points of tangency 
to the cutter surface during manufacture. The conditions 
necessary for envelope existence are given kinematically by the 
equation of meshing. This equation can be stated as follows: the 
normal of the generating surface must be perpendicular to the 
relative velocity between the cutter and the gear-tooth surface 
at the point in question (ref. 1). 

The coordinate transformation procedure that will now be 
described is required to locate any point from the head cutter 
into a coordinate system rigidly attached to the gear being 
manufactured. Homogeneous coordinates are used to allow 
rotation and translation of vectors simply by multiplying the 
matrix transformations. The method used for the coordinate 
transformation can be found in references 1, 5, and 8 to 11. 

Let us begin with the head-cutter coordinate system S ( shown 
in figure 3. This report assumes that the cutters are straight sided 
(not curved as commonly used on the wheel for final grinding). 
Surface coordinates u and 6 determine the location of a current 
point on the cutter surface as well as the orientation of the current 
point with respect to coordinate system S c . Angles \p ih and \ p ()h 
are the inside and outside blade angles. The inside and outside 
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Figure 3.— Head-cutter cone surfaces. 


blades cut the convex and concave sides of the gear teeth, 
respectively. A point on the cutter blade surface is determined 
by the following: 


r cot i p — u cos \J/ 


1 

where fixed value r is the radius of the blade at x c = 0, and \p 
is the blade angle. Parameters u and 6 locate a point in system 
S c and are unknowns whose value will be determined. 

The head-cutter coordinate system S c . is rigidly connected to 
coordinate system S s (fig. 4). System S s is rigidly connected to 
the cradle that rotates about the x m axis of the machine 
coordinate system S nr Coordinate system S m is a fixed 
coordinate system and is connected to the machine frame. 

To reference the head cutter in coordinate system S s , the 

following transformation is necessary: 
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y s y m 



Left-hand member; s = O c O s 



Figure 4.— Orientation of cutter, cradle, and fixed coordinate systems S c , S s , 
and S nr respectively. 


M = 

1Y1 n\s 


1 


0 


0 


0 


0 cos (j) c ± sin (j> c 0 

0 4= sin 4> c cos 0 


0 


0 


0 


1 


( 3 ) 


Coordinate system S m locates the machine center, and 
coordinate system S p orients the pitch apex of the gear being 
manufactured. The transformation from coordinate system S m 
to coordinate system S p requires the machine tool settings L m 
and E m along with dedendum angle <5 from the component 
design (see figs. 5 and 6). Machine tool settings L m and E m can 
be found from the summary sheet that typically accompanies a 
gear, or the methods in reference 8 can be used. Reference 8 
converts standard machine tool settings for the sliding base, the 
offset, and the machine center-to-back into settings L m and E m 
as shown in table I. The transformation matrix is given by 



M..,' — 


1 0 

0 cos q 
0 ±sin q cos q 
0 0 0 


0 0 
=Fsin q T s sin q 


( 2 ) 


^ cos q 

1 


where q is the cradle angle and 5 is the di stance between the 
Coordinate system S s and S c origins (s = O c O s ). The upper 
and lower signs preceding the various terms in this matrix 
transformation (and the rest of the paper) pertain to left- and right- 
hand gears respectively. 

Now, to transform from S s to the fixed coordinate system S m , 
the roll angle of the cradle <t> c is used. This transformation is 
given by 



Figure 5.— Coordinate system orientation to generate a right-hand gear surface 
= 0 shown here). 
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Front view 


Figure 6. — Coordinate system orientation to generate a left-hand gear surface 
(0 ( . = 0 shown here). 


COS 6 0 

0 1 


! Mp m 


-sin <5 

0 


sin <5 0 cos <5 

0 0 0 


-L,„sin 6 

±E m 

L m cos <5 


(4) 


This is shown in figure 5 for a right-hand member and in figure 6 
for a left-hand member. Figure 7 is given to clarify the orientation 
of the coordinate systems and machine tool settings (L m ,E m ). 

The next transformation involves rotation of system S p to S a . 
The common origin for coordinate system S p and S a locates the 
apex of the gear under consideration with respect to coordinate 
system S m . This requires rotation about y a by the pitch angle 
\jl (see fig. 7). This is given by 



cos n 0 sin /x 0 

0 10 0 
-sin /x 0 cos /x 0 

0 0 0 1 


(5) 


The final transformation is from coordinate system S cl to 
coordinate system S w , which is fixed to the component being 
manufactured. A rotation about the z^-axis through an angle </> w 
is required. Angle 0 H „ shown in figure 8, is the workpiece 
rotation angle; it is directly related to the angle of rotation of 
the cradle <t> c (this relationship will be described in the next 
section). The S a to S w coordinate transformation is given by 


TABLE I.— SIGN CONVENTIONS OF MACHINE-TOOL SETTINGS 
[From ref. 8.] 


Setting 

Sign 

Right-hand member 

Left-hand member 

Cradle angle, q 

+ 

Counterclockwise (CCW) 

Clockwise (CW) 


- 

Clockwise (CW) 

Counterclockwise (CCW) 

Machining offset, E tn 

+ 

Above machine center 

Below machine center 


- 

Below machine center 

Above machine center 

Machine center-to-back. 

+ 

Work withdrawal 

Work withdrawal 

Xmcb 

- 

Work advance 

Work advance 

Sliding base, X SB 

+ 

Work withdrawal 

Work withdrawal 


- 

Work advance 

Work advance 

Vector sum of X SB 

+ 

XSB* + * XmcB' ~ 

^SB’ +’ XmCB' ~ 

and X MCB 

— 

%SB’ ~ XmCB’ + 

X$B' ~ *' X M CB' + 
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• Plane n 



z a> z w 


(b) 


(a) Machine settings and orientation. 

(b) Plane ir and orientation of generated gear coordinates. 


Figure 7.— Orientation of machine settings and generated gear coordinate systems. 





Right-hand member 

Figure 8.— Rotation of workpiece for left- and right-hand gears during tooth- 
surface generation. 
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COS 0 M . 

±sin 4> w 

0 
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=Fsin 4> w 

cos </>„. 

0 

0 

(6) 

Tooth Surface Coordinate Solution Procedure 

0 

0 

1 

0 


In order to solve for the coordinates of a spiral bevel gear- 

0 

0 

0 
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tooth surface, the following items must be used simultaneously: 


Using these matrix transformations, we can determine the 
coordinates in S w of a point on the generating surface from 


r„. = [AU [M ap \ [MpJ [M ms ] [M S( ] r ( 


(7) 


or 


or 


the transformation process, the equation of meshing, and the 
basic gear design information. The transformation process 
described previously is used to determine the location of a point 
on the head cutter in coordinate system S w . Since there are 
three unknown quantities («, 6 , and 0<), three equations 
relating them must be developed. 

Values for u, 6 , and </>,. are used to satisfy the equation of 
meshing given by references 1 and 9: 


r M . = [MJMc) ] 1^] [M pm ] [M m (<l> c )][M sc ]r c (u,0) (8) 


n • V = 0 


(10) 


r tt . = r Au,0,(t> c ) 


(9) 


This transformation describes the location of a point in the 
gear fixed coordinate system based on machine settings (L m , 
E nv q, s, r y and yp), parameters (u, 6 , and 0 r ), and gear 
design information (fjL and 5). 


where n is the normal vector to the cutter and workpiece 
surfaces at the specified location of interest, and V is the 
relative velocity between the cutter and workpiece surfaces 
at the specified location. From the reference 9 equation of 
meshing for straight-sided cutters with a constant ratio of roll 
between the cutter and workpiece, equation (10) is defined as 
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(w — r cot 0 cos i/Ocos 7 sin r 


+ s[(m cu . — sin 7)cos^, sin 0 =F cos 7 sin \p sin (<7 — 0 ( )] 

± £ w (cos 7 sin 0 + sin 7 cos \p cos r) 

- L,„ sin 7 cos 0 sin r = 0 (11) 

where 7 is the root angle of the component being 
manufactured, and 


t = (6 T q ± </> ( .) 


and 


m n . 


AO 


>•) 


( 12 ) 


(13) 


where m cw is the ratio of angular velocity of the cradle to that 
of the workpiece. Since the ratio of roll in this report is 
assumed to be constant, equation (13) can be written as 


x a , x w 



View A 



Figure 9.— Orientation of gear to be generated, with assumed positions / and 


00 


(u) 


m nv 


and 


The axial position must match the value found from 
transforming the cutter coordinates S c to workpiece 
coordinates S w . This is satisfied by the following (fig. 9): 


a {w) =— w (f) = — 
dt dt 

therefore 


<t>c 

d> w = — (14) 

m cw 

Equation (14) is the relationship between the cradle and 
workpiece for a constant ratio of roll and is used directly in 
equation (6). 

Gear design information is then used to establish an 
allowable range of values of the radial (r) and axial ( z ) 
positions that are known to exist on the gear being generated. 
This is shown in figure 9. 

First the equation of meshing must be satisfied. This was 
shown earlier to be 

n • V = 0 
or 

f\(u, 6, 0,) =0 (15) 


z w -z = 0 (16) 

or 

f 2 (u, 0, 4> c ) = 0 (17) 

Finally the radial location from the work axis of rotation must 
be satisfied. This is accomplished by using the magnitude of 
the location in question in the x w -y w plane (see fig. 9): 

r - [xl + yl)°' 5 = 0 (18) 

or 

Mu, e, (t> c ) =0 (19) 

Now a system of three equations (eqs. (15), (17), and (19)) 
is solved simultaneously for the three parameters u, 0, and 
0 r , for a given gear design with a set of machine tool 
settings. These are nonlinear algebraic equations that can be 
solved numerically with commercially available mathematical 
subroutines. These equations are then solved simultaneously 
for each location of interest along the tooth flank, as shown 
in figure 10. In the SURFACE program a 10 by 10 grid of 
points is used on each side of the tooth. From the surface grids, 
the active profile (working depth) occupied by a single tooth 
is defined. 
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Figure 10.— Calculation points (10 by 10 grids, i.e., 100 points each side) 
for concave and convex sides of tooth surface. 


Application of Solution Technique 

An application of the techniques previously discussed will 
now be presented. The component to be modeled was from 
the NASA Lewis Spiral Bevel Gear Test Facility. A photo- 
graph of the spiral bevel gear mesh is shown in figure 1 1 , and 
the design data for the pinion member are shown in table II. 

The gear design data were used along with the methods of 
reference 9 to determine the machine tool settings for straight- 
sided cutters (see table II). These values were then used as 
input to SURFACE. This program calculates the coordinates 




Figure 1 1.— NASA spiral bevel gear test rig components. 


of the concave and convex sides of the gear tooth (fig. 10), 
orients the surfaces such that the top land is of the proper 
width, and then generates the required data for the three- 


TABLE II.— EXAMPLE CASE OF SURFACE COORDINATE GENERATION 
(a) Pinion design data 


Number of teeth 

pinion 12 

gear 36 

Dedendum angle, deg 1.0 

Addendum angle, deg 3.883 

Pitch angle, deg 18.433 

Shaft angle, deg 90.0 

Mean spiral angle, deg 35.0 

Face width, mm (in.) 25.4 (1.0) 

Mean cone distance, mm (in.) 81.05 (3.191) 

Inside radius of gear blank, mm (in.) 15.3 (0.6094) 

Top land thickness, mm (in.) 2.032 (0.080) 

Clearance, mm (in.) 0.762 (0.030) 


(b) Generation machine settings 



Concave 

Convex 

Radius of cutter, r, mm (in.) 

75.222 (2.9615) 

78.1329 (3.0761) 

Blade angle, 0, deg 

161.358 

24.932 

Vector sum, L,„, mm (in.) 

1.0363 (0.0408) 

-1.4249 (-0.0561) 

Machine offset, E nv mm (in.) 

3.9802 (0.1567) 

-4.4856 (-0.1766) 

Cradle to cutter distance, s , mm (in.) 

74.839 (2.9646) 

71.247 (2.8050) 

Cradle angle, q , deg 

64.01 

53.82 

Ratio of roll, m cw 

0.308462 

0.321767 

Initial cutter length, a , mm (in.) 

239.5 (9.43) 

181.1 (7.13) 

Initial cutter orientation, 6, deg 

120.0 

120.0 

Initial cradle orientation, 0 t ., deg 

0.0 

0.0 
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dimensional modeling program PATRAN (ref. 31). The details 
of the procedure are described in the following paragraphs. 

Surface Coordinate Calculation 

Using figures 10 and 12 as references, we will describe the 
calculation procedure for surface coordinates. First, the 
concave side of the tooth is completely defined before moving 
to the convex side. These points are calculated by starting at 
the toe end and at the lowest point of active profile. Nine steps 
of equal distance are used from the beginning of the active 
profile to the face angle (addendum) of the gear tooth, and 
then back to the next axial position (see fig. 12). The procedure 
is repeated until the concave side is completely described. Then 
the same procedure is followed for the convex side. 

In the discussion of surface coordinate determination, the 
cutting blades were described as straight sided. The point 
radius r (eq. (1)) is the radius the cutter would have if it were 
projected down to the y c -z c plane (figs. 3, 5, and 6). Actually, 
the blade has a theoretical point width and corner radii that 
generate the portion of the tooth from the working depth to 
the root cone (see fig. 13). In the section Coordinate Transfor- 
mations, this part of the cutting blade was not modeled, so 
the current analysis by SURFACE either assumes a full fillet 
radius between the lowest point of active profile on adjacent 
teeth or sets the fillet radius equal to the clearance (fig. 14). 

Concave and Convex Orientation 

Since both sides of the gear-tooth surface are not analyzed 
simultaneously, for proper alignment of the surfaces, their 
orientation relative to each other must be established by 
determining the amount of rotation in the fixed coordinate 



Figure 12.— Cross section of calculation grid: I\ root angle: 6, dedendum 
angle: addendum angle: r jt inside radius of gear blank. 


system S w . This is done by checking the tooth thickness at the 
face angle on the toe end of the gear tooth. (Also, in the case 
of a gear, remember that a given cutting operation using cutters 
as shown in fig. 3 actually cuts adjacent teeth on the convex 
and concave sides simultaneously.) The distance between these 
two locations must correspond to the top land width. The 
convex surface is then rotated according to the angle deter- 
mined by the points at the face angle at the toe position. This 
is shown in figure 15. (Note that this same procedure could 
have been done by considering the tooth mean circular 
thickness instead of the tooth thickness at the face angle as 
was done in this report.) 


*c 



Figure 13.— Detailed view of straight-sided cutter geometry. 
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Working 



Figure 14.— Two types of fillet and root radius regions used by SURFACE 
program; p c , constant fillet and root radius; pp fillet radius equal to the 
clearance; N , number of teeth; r h inside radius of gear blank. 


Yw 



Figure 15.— Orientation of concave and convex sides of gear-tooth surfaces 
to attain proper top level width by rotation of convex side of tooth; P,, 
concave side location of face angle point at toe end of tooth; P 2 and P 2 , 
initial and final convex side locations, respectively, of face angle point at 
toe end of tooth; £, angle through which P 2 is rotated; and w (l and w h , 
desired and initial top level thicknesses, respectively. 




Figure 16.— Hidden line plot of one-tooth model generated by the SURFACE program using PATRAN (constant fillet and root radius model). 
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Generation of Three-Dimensional Model Data 

Once the surface coordinates are described and properly 
oriented, then the data necessary for PATRAN have been 
produced. Currently, the model in SURFACE can have two 
different Fillet and root radii configurations. The fillet radius 
can be constant at the cross section between adjacent teeth or 
be equal to the clearance with a flat between the fillet radii 
of the adjacent teeth at the root angle. Also a constant inside 
radius of the gear blank is used in this modeling method. These 
assumptions are depicted in figure 14. 

At this point, we have produced a one-tooth model for use 
in PATRAN. Now, the analyst must determine how complex 
the model need be for a given application. If a complete gear 
is required, simply rotate the one-tooth model in PATRAN. 

Once the required number of teeth have been described, then 
the finite element mesh density and the boundary condition 
information are generated within the PATRAN environment. 
PATRAN produces the bulk data deck for MSC/NASTRAN 
and many other computer codes. The example given in this 
report used MSC/NASTRAN for the static analysis. 


Example Model and Results 

From the one-tooth model described earlier the analysis 
techniques can be demonstrated. The model shown in figure 
16 is that for a constant fillet and root radius (also called full 
fillet) model. The fillet and root radius on the convex side has 
been added along with the tooth section (without the tooth) 
to make the model symmetric about the tooth centerline. Figure 
16 shows a hidden line plot of the finite element mesh with 
eight-noded isoparametric three-dimensional solid continuum 
elements. This model has 765 elements and 1 120 nodes. The 
boundary conditions are shown in figure 17. A 1724-MPa 
(250 ksi) constant pressure load was applied normal to the 
tooth surface of nine elements, and the two edge surfaces of 
the gear rim had all degrees of freedom constrained. 

The results were calculated by MSC/NASTRAN and were 
subsequently displayed by PATRAN. Figure 18 shows the 
principle stresses, and figure 19 shows the total displacements 
for the boundary conditions shown in figure 17. 


1724 MPa (250 ksi) 
pressure load 



Figure 17. — Boundary conditions for constant fillet and root radius model for the example application. 
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Stress, 

ksi MPa 



Figure 18.— Major principal stress for boundary conditions specified in figure 17. 


Displacement, 
in. pm 

4.08x1 0" 3 



Figure 19.— Total displacement for the boundary conditions specified in figure 17. 
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Summary of Results 

A method has been presented that uses differential geometry 
techniques to calculate the surface coordinates of face-milled 
spiral bevel gear teeth. The coordinates must be solved for 
numerically by a simultaneous solution of nonlinear algebraic 
equations. These equations relate the kinematics of manu- 
facture to the gear design parameters. Coordinates for a grid 
of points are determined for both the concave and convex sides 
of the gear tooth. These coordinates are then combined to form 
the enclosed surface of one gear tooth. A computer program, 
SURFACE, was developed to solve for the gear-tooth surface 
coordinates and provide input to a three-dimensional geometric 
modeling program (i.e., PATRAN). This enables an analysis 
by the finite element method. An example of the technique 
was presented. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, December 19, 1990 
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